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1 Introduction 

The discrete nonlinear Schrodinger equation (d-NLS) 

Aix + h,t) + Aix -h,t) -2A(x,t) ,,,2, ^ , 
lAt +p— ^-^ — ^-^+q\A\^A = 0,pq ^0, p and q real(l.l 

occurs in many physical disciplines, but unfortunately is not integrable. Hence most inves- 
tigations must be conducted via numerical computations or other means of approximation. 
In contrast, the Ablowitz-Ladik (AL) model [1], 

^A, + + ^) + - + + ^' ; - ^' \Aix, t)P = 0,(1.2) 
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possesses A^-soliton solutions, periodic waves, an infinite number of conservation laws and 
a host of other properties associated with integrable systems [271 [2H]- Although d-NLS 
occurs much more frequently than AL in physical situations, there are still some situations 
where AL is relevant, e.g., higher order evolution equations [T7], optics [23) . and protein 
dynamics [21] . Indeed generalizations of AL systems have received quite intensive attention 
in the literature [24j . 

The goal here is to investigate a class of coupled (or vector) Ablowitz-Ladik systems 
with saturable nonlinearities. By "saturable" NLS-type equation, we mean that, in the 
nonlinear term, the coefficient of A, typically \A\'^, is replaced by some expression whose 
limit is finite when |^| — > +oo, such as |A|^/(1 + \B\'^) or, in the case of two-component 
NLS, byl^p/d^p + lAp). A brief and certainly incomplete review of some existing works 
in the literature will provide the motivation of the present work. 

We first look at recent advances on the AL system. Straightforward AL extensions of 
vector nonlinear Schrodinger equations exhibit solitons with elastic collision properties [2] . 
Indeed collision properties of chains of near identital solitons have been studied too [3 [9]. 
Other variants, e.g. coupled discrete evolution systems with nonlinearities resembling the 
AL model [4i|22j, have been considered. However, they are usually not integrable, and 
propagation of solitons may be unstable. 

A remark on saturable nonlinearity is also in order. Saturable nonlinearities probably 
first arose in the setting of nonlinear optics. In an optical fiber, the Kerr nonlinearity 
typically generates the conventional nonlinear Schrodinger equation (NLS). However, for 
short pulses and a high input peak pulse power, these Kerr nonlinearities cannot ade- 
quately describe the field-induced change in the refractive index. In such a case, it is 
sufficient to introduce a saturation term in the nonlinearity, leading in the simplest case 
to the one-component saturable (continuous) nonlinear Schrodinger equation [15| 

iAt + pA^^ + q ' = 0, (1.3) 

i + /ii^r 

in which p, q, fj, are real constants. Other types of nonlinearities have also been considered 

m. 



Its d-NLS-type discretization 

A(x + h,t) + A(x-h,t)-2A(x,t) lAl'^A , , 
^ ' ^ + q- ] = 0,pqu + 0, 1.4 

in which the saturation coefficient v is dimensionless, has a considerable importance in 
nonlinear optics [11], but very few analytic results are known. Apart from plane waves and 
some perturbative solution [26], the main particular solution so far is an elliptic stationary 
wave solution [19] , recalled as Eq. (j3.5|) in section [3l In particular the problem remains 
open to find an exact travelling wave of ()1.4p with an arbitrary velocity. 

When the saturation term 1/(1 + z/(g/i^/p)|yl|^) is replaced by an arbitrary analytic 
function of \A^ and, at the same time, the nonlinearity is partly discretized in the manner 
of Ablowitz and Ladik, 

A{x + h,t) + A{x -h,t) - 2A{x, t) 



iAt + P- 



/l2 



+ ,(i.(|^|)ii(^±M±jl(^J^ + G(|A|)A(x,t)) =0, (1.5) 
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one can require the existence of a special class of moving kinks and pulses and thus 
determine F and G [13j, with the main result that Eq. (jl.4p admits no stationary pulse 
in the considered class. A similar inverse approach has been followed in Ref. jl8j . 

Over the years, dynamics of pulses in saturable media has continued to attract at- 
tention. An example of a coupled waveguide, which we plan to expand on here, focuses 
on a saturable, nonlinear medium with cross phase modulation. The propagation of two- 
component spatial optical solitons carrying angular momentum has been investigated |25j . 
The reduction of this system to one transverse dimension defines the two-component sat- 
urable NLS 



^ . \Af + ai\B\^ 

lAt + AiMxx + q . A = 0, 

1 + (\A\^ + \B\^ 

iBt + hpBxx + q 7 z^B = 0, 

1 + S2 1^1 



(1.6) 



in which p, q, Xi, X2, si, S2, cri, 02 are real constants, and this system has been shown [6] to 
admit at least one stationary solution (of elliptic type) with a nonconstant total intensity 
-|- l-Bp, so that the system is truly saturable. In section [5l we extrapolate this solution 
to a travelling wave having an arbitrary velocity. Another two-component saturable NLS, 



iAt + Aip^xx + q 7 — TV = 0' 

1 + si {\A\^ + \B\^ 



B 

iBt + X2pB^x + q J = 0, 

1 + S2 1^1 



(1-7) 



describes optically induced nonlinear photonic lattices in which discrete solitons have been 
observed pTTIfTi]. 

From the perspective of nonlinear science, as well as from the viewpoint of applications, 
it is therefore instructive to consider NLS-type equations which are discrete in space and 
possess a saturable nonlinearity. Given the elegant nature of the AL system, it seems 
plausible to study such a discrete saturable NLS with the nonlinear term discretized in 
the manner of AL rather than of the d-NLS. In fact preliminary studies of such systems 
have already begun |T9| [18]. 

The structure of this paper is the following. In section [2] we explain how to perform 
computations which minimize the algebraic complexity. In section [3l we first describe a 
new family of stationary solutions for a single component AL with saturable nonlinearity. 
These new solutions are different from those given in the literature |19l llSj . They are 
delineated by a cyclic combination of elliptic functions in the compact and efficient notation 
due to Halphen. 

Then in section [U we consider a system of two coupled saturable continuous NLS and 
generalize a structure previously obtained [6j to an arbitrary velocity. 

Finally in section [5] we introduce an AL-type discretization of system (II. 6p and find 
exact solutions with a nonconstant total intensity \A\'^ + \B\^, so that the satur ability 
feature is effective. These solutions are stationary and of the elliptic type. To help the 
reader, we have gathered in Appendix all the mentioned solutions, whether existing or 
new. 
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2 Practical computation 

Since the twelve Jacobi elhptic functions [3l Chap. 16] are equivalent in the complex plane, 
any elliptic solution can be presented as twelve equivalent complex expressions. However, 
in physics we are only interested in those expressions which are bounded on the real axis, 
i.e. which only involve the six functions sn, cn, dn and sd, cd, nd of the canonical argument 
{Qx,k), in which the Jacobi modulus k lies between and 1 and the Jacobi nome Q is 
positive P §18.9.9, 18.9.12], 

m = k'^ = ^ q2 ^ ei - 63, < m < 1, < Q^ eg < 62 < ei. (2.1) 

ei - 63 

In order to minimize time and effort, we will follow a two-step procedure. 

1. Firstly, to find a complex (mathematical) solution A = f{haix),h^{x),h.'y{x)), in 
which (a, (5, 7) is an arbitrary permutation of (1, 2, 3) and (hi, h2, hs) are the elliptic 
functions introduced by Halphen |16j . defined by the differential system 



^ = -h/3h„h^ - 



h/3 h7, hl-hl = -ep + e^, ei + 62 + 63 = 0, (2.2) 



with the choice of sign \mix^QX\ia{x) = +1. The invariance under any permutation 
of (a, /?, 7) is an important practical advantage of the Halphen notation over the 
Jacobi notation. 

2. Secondly, for each such complex solution, to convert the Halphen trio (hi,h2)h3) to 
both bounded Jacobi trios (sn, cn, dn) and (sd, cd, nd) by formulae such as 

hi (re - xo) = ai dn{Qx, k), h2{x - xq) = ^2 cn{Qx, k), 
hsix - xo) = assn{Qx,k), (-23) 

222 "203 ^,2 

fli = 63 - ei, a2 = 63 - 62, 03 = 62 - 63, = Qk , 

«i 

in which xq is an immaterial complex shift. Full details can be found in Ref. [5l 
App. B]. 

3 One-component saturable discrete NLS 

For convenience, let us introduce the constant ao, a kind of unit of amplitude, 

al = -2^. (3.1) 

Let us first consider Eq. (jl.4p . 

A(x + h,t)+A(x-h,t)-2A(x,t) lA^A ^ , ^ rr^ 

l + v{qh^/p)\A\ 

If one assumes for a solution of (jl.4p the expression (first step as defined in section [2]) 



A = ai 



, , h/3(x-ei) 

ha(a; - ct) + lOi— — 

h7(x - ct) 



^i{K{x-ct)-u}t) ^ (•3_2) 
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one finds six constraints among the seven introduced constants (ai, 61, c, to, K, e^, C/j), 



c = 0, CO = I 2 , tan Kh 

bi{e-y - e/3){e^ - - bj) = 0, 



2 2 

ai = a, 



2v 



' °4z.2/,2 cos Kh{hl{ep - e^) + h^(/i)) ' 

h^(/i) [6?(h2 (/i) h^(/i) - (e^ - e«) + hi{h) h6(/i) 



(3.3) 



h;3(/i) cos Kh {bl{ep - e^) + h^(/i))' 



The factorized form of the fourth constraint (quite similar to that obtained in [7]) defines 
three solutions, each depending on one arbitrary constant and requiring that v be nonzero. 
In order to write them explicitly, let us introduce the nonzero constant 



= h^(/i) - (e„ - e^)(e„ - e^) = 2ha(/i) hpih) h^{h) ha(2/i). 
These three solutions are 
o [hi = 0): the elliptic solution already obtained |19j . 



(3.4) 



A = aiha(x)e ai = a, 



2\ip{h)h^{h) 



2v 



h/3(/l) h^(/i)^ 



(3.5) 



o (e^ — = 0): a dark one-soliton solution 
A = ai 



tani^/i = -6i-tanh — 
2 



— tanh —X + zoi 
2 

•^1 ~ ^0 ( ~~r I ^i'^^ — cosh — COS A /i, 2z^ = cosh — cos Kh^ 
y /^/i y 2 2 2 



= 0): a new elhptic solution, diff'erent from the previous one, 



(3.6) 



A = ax 



ha(x) + ih 



h/3(a 



,2 _ „2h/3(^)h7(/i)A^7 



„i(/^x-c^i)^ ^2 ^ 
h7(/l)A^a 



(3.7) 



h'^N'^cosKh 



hl3{h)N^ cos ET/i 



In the continuum limit, the coefficient goes to 1/2 for all three solutions, so the limit 
for Lo requires expanding v to second order in h. Each solution goes to some solution of 
the continuous NLS, 

o (61 = 0): 



6/3, e-y = arbitrary, A = aoho(2;)e , 2^ = 1, ui = —3pea, 



(3.8) 
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(e^ - e/3 = 0): 



bi,K = arbitrary, A = 



K K 

— tann —x + zoi 



2u = l, L0= {3bi + —]p, 



^i{—bix—ujt) 



(3.9) 



(e^ -ea-bl = 0): 



6/3, = arbitrary, ^ = ag 
2z^ = 1, u! = (4e/3 + 5e^)^j. 



ha(x) + ibi 



h/3(a:) 



(3.10) 



Following the second step of section [21 various bounded solutions represented by the 
complex expression (13. TP can be found in terms of bounded Jacobi functions, such as 



A = (sn +i cd)e 



i{Kx—ujt) 



, (cn+zsd)e 



i(Kx—u)t) 



(3.11) 



we leave this exercise to the interested reader. 

On this NLS limit, it is easier to give the precise difference between the two elliptic so- 
lutions (|3.8p and (|3.10p . i.e. between (|3.5p and (|3.7p . Indeed, the travelling wave reduction 
of NLS 



A{x, t) = yM(6e^^"^* + C = x-ct, 

admits as general solution the elliptic functions 

M= =ag(p(^- ^0,52,53) -eo), ^' 



(3.12) 



2^ + "° 



iijp - 2 



(3.13) 



-^(4eo -5260 -53), 



therefore the above two elliptic solutions can only be particular cases of it. The first 
solution (j3.8p is the well known Lame function corresponding to cq = ea,Ki = 0. As to 
the second solution ()3.10p . by establishing the ODEs for |Ap and (avgA)' one finds the 
values 



eo 



3£'q, + Ejs, Ki — —2aQyJ Ej — Ea{Ea — E, 



(3.14) 



in which the elliptic invariants Ej of (j3.10p are related to those Cj of (j3.8p by the Landen 
transformation which divides by two only one of the two periods |12| §13.23], 



p{z\u,uj') -> p(^z 



iO 



(3.15) 



This is the nonzero value of Ki which proves that both solutions are indeed different. 

We have checked that the saturable discrete NLS (|1.4p does not admit a solution de- 
pending on two elliptic functions with different periods such as the one [5] admitted by 
the discrete integrable NLS of Ablowitz and Ladik. 
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4 Two-component saturable continuous NLS 

We next consider the two-component saturable continuous NLS Eq. (jl.6p . 



iAt + XipA^^ + g- 



l + si( \A\'^ + \B\ 



iBt + \2pBxx + q 



^ = 0, 



-8 = 0, 



In the continuum case, when one requires the total intensity \A\'^ + \B\'^ to be nonconstant, 
so that the system (|1.6|) truly display saturability, then only one elliptic solution is known 
[6], in which the moduli |^|, \B\ have a zero velocity c. However, if one assumes 



A = ai hp{x - ct) h^{x - ct)e^(^«(^ " " 

B = a2 hpix - ct) ha{x - ct)e'(^bix - ct) - iv^t) ^ 



(4.1) 



the above mentioned solution, which contained no arbitrary parameter at all, is extended 
to a new solution with one arbitrary parameter (the velocity c), defined by the algebraic 
relations 



c = arbitrary, of = a, 



+ 



-a, 
A2 



a = 9a| 

C'-y I 

Ka 



1 — fJi 1 — (72 



1 



1 



1 / 1 — fJl 1 — (72 

Aisi A2S2 2 



Aisi 

1 - (72 



1-ai 



SagAiSi 3aoA2S2 



A2S2 



2pXi 
c 



1 ^2 „2x2 



1-^ 

9a^S2 



2 \2 



+ 2- 



a^A 



(4.2) 



UJb 



Ka+ 2 



2pA2' 

SpAiCc 



2P o„x „ ISpAfag 



a(l - (7i) 
ISpAlag 



-KbH 2-3pA2e^+ . 

2 52^0 "(1 - ^2) 



and there exists exactly one constraint among the fixed coefficients, 

0, 



1 - CJl _^ 1 - fT2 



AiSl A2S2 

which implies si ^ S2- In particular, the dispersion relations take the form 



-Kj + constant,-, j = a, b. 
2 



(4.3) 



(4.4) 



Since the three roots Ca are always real, this complex solution defines three bounded 
solutions (in terms of sn, cn, dn), identical to those mentioned in [6, Eqs. (12)-(13)], 
i.e. {A,B) = (sndn, cndn), (cnsn, dnsn), (dncn, sncn), with minor adjustments for the 
additional dependence on c. 
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The long wave hmit of this solution is {e/3 — ea)[ep — e^) = 0, i.e. 

J ^ = aisechA;(x - ct)tanhA;(x - ct)e*'^^«(^ ~ "^^^ - ^J) 
\ B = a2 sech2 k{x - ct)e'(^b{x - ct) - uj^t) ^ 

with the relations 

c = arbitrary, of = = 



A2S2 



3%S2 \Xi A2 



Ka 



2pX, 



2pX2 

— 2 ( 1 - — I ' ^fe 



c ip 



4pXi 6s2% V A2 
and the two fixed constraints 



1-0-1 = — 1-— , 1-0-2 = 1- —. 

S2 \ X2 Ai 



+ 



4pA2 3s2aQ V Ai 



A2 

1 - -A 



(4.5) 



(4.6) 



(4.7) 



5 Two-component saturable discrete NLS of A-L type 

Let us finally consider an AL-type discretization of system (jl.6p , namely 

A(x + h,t)+ A{x -h,t) - 2A{x) 



iAt + Xip- 



+ q 



\A\' + ai\Bf^ 



A{x + h,t) + A{x - h,t) 



l + s,[\A\^ + \B\^ 



0, 



iBt + A2P 
+ q- 



B{x + h,t)+ B{x -h,t)- 2B{x) 

\B\^ + a2\A\^ B{x + h,t)+B{x-h,t) 



(5.1) 



l + S2{\A\^ + \B\^ 



The expression (j4.ip where x is now a discrete variable x = nh is still a solution of (j5.1 
but with much less freedom among the parameters. It is defined by 



c = 0, of = a, 02 = —a, Ka = 0, Ki, = 0, 



si = + 



S2 



a{ea - e^)(h^(/i) - Se^) 
1 

a(e^ - ea)ih.aih) — 3e^) 



, Aisi 
, A2S2 



al hUh){3liUh) +ea-ei3)' 



L0„ 



Wfe = + 



2pa{e^ - e„) \ip{h){^{h) - 3e„)(h|(/i)(h/3(/i) - h^(/i)) + 2(e„ - ep) \i^f{h)) 

al hl{h){?,\il{h) + ea - ep) ' 
2pa{ea - e^) hf3ih){hl{h) - 3e^)(h|(/i)(h/3(/i) - ha(/i)) + 2(e^ - ep) ha{h)) 



(5.2) 



^1 = . 9 /7 N/o, 9 /, \ , i — 0-2 — 
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The velocity c now vanishes (while it is arbitrary in the continuous case) and the fixed pa- 
rameters Ai, A2, (Ji, (72, si, S2 obey three constraints (as compared to one in the continuous 
case). In particular the relation 



1 

si 



1 

S2 



-2a(ec 



^7/ ) 



(5.3) 



requires si and S2 to be different if the elliptic functions remain nondegenerate. 

In order to remove the constraint c = 0, we have tried the assumption (inspired by the 
relations (lOI)! 



A 
B 



-Ol 

-a2 



h.(g + V2)-h.(g-V2) j(K„£ - u^^t) ^^ = ^_^t, 
h^(g + h/2) - h^(g - h/2) - u.t) ^ ,^ 



(5.4) 



but the system admits no such solution at all. 

The long wave limit of the solution (j4.ip . (j5.2p is (e/3 — ea){ef3 — e-y) = 0, and it leads 
to the same expression (14. 5p as in the continuous case. 



6 Discussion and conclusion 

A system of coupled discrete evolution equations, with cubic nonlinear terms of the 
Ablowitz-Ladik type and the presence of saturable nonlinearities, has been studied. For 
the single component case, a new type of periodic waves is formulated in terms of the 
elliptic functions introduced by Halphen, which maintain a high degree of symmetry. This 
family of periodic patterns is different from those established earlier in the literature. 

For the two-component case, solutions in terms of single and products of elliptic func- 
tions are found, which both correspond roughly to sinusoidal structures. As expected in 
any complicated nonlinear system, significant constraints on the properties of the sys- 
tem, i.e. coefficients of dispersion, cross phase modulation, and saturation, must exist for 
analytical progress to be made. 

A very important issue which has not been addressed in the present work is the stability 
of these wave patterns. Computational tests must be conducted which will be performed 
in the near future. 
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Appendix. Summary of solutions 

Table 1. Various types of NLS equations considered in the text. Types= (discrete or continuous), 
(one- or two-component), (saturable or not), (AL or non-AL). Y=yes, N=no, B=both, ell=elhptic, 
tri=trigononietric. A star (*) denotes apparently new results. 
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